Aparici Seguer atal.\ Moiuscos infralitorales de Pinedo-El Saler, Valencia 



Estaciones 


Figura 3. Representacion grafica de los valores de diversidad para cada estacion. A: fndice de 
Riqueza Especffica No; B: fndice de Diversidad H'; C: fndice de Homogeneidad E5. 

Figure 3. Graphic representations of the diversity values in each station. A: No Richness index; 
B: H' Diversity index; C: E5 Evenness Index. 


1 0,75 0,50 0,25 0 



Figura 4. Dendrograma de afinidad entre estaciones, construido segun el algoritmo UPGMA. 
Figure 4. Dendrogram of stations classification, constructed according to the UPGMA method. 
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desde especies caracteristicas exclusivas 
de la Comunidad de Arenas Finas Su- 
perficiales (SFS), como Donax semistria- 
tus, Tellina tenuis, Lentidium mediterra- 
neum y Chamelea gallina, hasta la Comu¬ 
nidad de Arenas Finas Bien Calibradas 
(SFBC), con Mactra stultorum, Tellina pul- 
chella, T nitida, Spisula subtruncata, Pan¬ 
dora inaequivalvis y Loripes lacteus. 

Gola es el unico transecto en que se ha 
registrado una disminucion de la salinidad 
en superficie (37,78%o), debida al impor- 
tante aporte de aguas continentales de 
drenaje de 1' Albufera. Este hecho no parece 
afectar a los moluscos alii estudiados. 
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Microdiferenciacion genetica de Mytilus edulis y 
M. galloprovincialis en el golfo de Vizcaya 

Genetic microdifferentiation of Mytilus edulis and M. 
galloprovincialis in the bay of Biscay 


Angel S. COMESANA y Andres SANJUAN* 


RESUMEN 

Mytilus edulis y M. galloprovincialis son dos formas de mejillones que habitan en las cos¬ 
tas europeas, coexistiendo e hibridando en proportion variable. Se estudiaron poblacio- 
nes de mejillones de la zona hibrida de ambos Mytilus en el Golfo de Vizcaya (Capbre- 
ton) en una zona expuesta al oleaje y en otra resguardada cercana, y en el nivel alto y 
en el bajo de fijacion. Los mejillones se analizaron electroforeticamente para cuatro 
genes enzimaticos diagnostics entre ambos Mytilus : esterasa-D, leucin aminopepti- 
dasa-1, manosa fosfato isomerasa y octopfn deshidrogenasa. 

Se encontro que los alelos caracteristicos de M . galloprovincialis para los cuatro genes 
son significativamente mas abundantes en las poblaciones de la zona expuesta de fija¬ 
cion (70 a 85%) que en las de la zona resguardada (20-70%); y en estas, mayor en las 
del nivel alto (50-70%) que en las del nivel bajo (20-25%). Se confirman y generalizan asf 
dos fenomenos anteriormente detectados en grado diverso en las Islas Britanicas. Se dis¬ 
cute la posible interaccion de los factores exposition al oleaje y nivel de fijacion, para 
comprender la composition genetica y la dinamica de las poblaciones mixtas de Mytilus. 
Ademas, la importante microdiferenciacion genetica en un area que bien pudiera ser la 
zona de contacto de M. edulis y M. galloprovincialis mas antigua y en la que cabria espe- 
rar una amplia homogeneizacion genetica, suscita la existencia de alguna ventaja adap- 
tativa todavia desconocida de M. edulis para poder mantenerse como tal en dicha area. 

ABSTRACT 

Mytilus edulis L. and M. galloprovincialis Lmk. are two forms of mussels which inhabit the 
European coasts, where hybridization between them occurs in variable proportions. 
Samples of mussels within the hybrid zone in the Bay of Biscay (Capbreton) from closely 
adjacent exposed and sheltered locations and, inside each one, from high and low sho¬ 
res were studied. Four enzyme genes, which are partially diagnostic between both Myti¬ 
lus, were analyzed by electrophoresis: esterase-D, leucine aminopeptidase-1, mannose 
phosphate isomerase and octopine dehydrogenase. 

The characteristic alleles of M. galloprovincialis for the four genes are significatively 
more abundant in the exposed populations (70 to 85%) than in sheltered ones (20 to 
70%), and in the latter, in high shore populations (50 to 70%) than in low ones (20 to 
25%). So these two phenomena, previously detected in variable degree in the British 
Isles, were confirmed and generalized. The possible interaccion between both factors, 
exposure to wave action and level of attachment, are discused to understand the genetic 
composition and the dynamics of the mixed populations of Mytilus. On the other hand, 
the important genetic microdifferentiation found in this area, which could be the oldest 
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contact area of M. edulis and M. galloprovincialis and where an extensive genetic homo¬ 
genization could be expected, suggests an unknown adaptive advantage of M. edulis 
allowing its persistance as a different form in this area. 


PALABRAS CLAVE: Mytilus edulis , Mytilus galloprovincialis , polimorfismos proteinicos, genes enzima- 
ticos, genetica de poblaciones, zona hibrida, diferenciacion genetica, Golfo de Vizcaya 
KEY WORDS: Mytilus edulis , Mytilus galloprovincialis, protein polymorphisms, enzyme genes, popula¬ 
tion genetics, hybrid zone, genetic differentiation. Bay of Biscay 


INTRODUCCION 

Mytilus edulis Linne, 1758 y Mytilus 
galloprovincialis Lamarck, 1819 son dos 
formas de mejillones que habitan las 
costas europeas y cuyas areas de distri- 
bucion en el litoral europeo surocciden- 
tal se ha determinado recientemente. M. 
galloprovincialis se encuentra desde el 
Mar Mediterraneo hasta el Canal de la 
Mancha y las Islas Britanicas, y M. edulis 
se extiende desde el litoral atlantico 
norte hasta tan al sur como la frontera 
hispanofrancesa (Fig. 1; SanjuAn, 
Quesada, Zapata y Alvarez, 1986, 
1990; Coustau, Renaud y Delay, 1991; 
SanjuAn, 1992; para revision ver 
Gosling, 1992a). Los estudios de genes 
enzimaticos y de polimorfismos de frag¬ 
ments de DNA mitocondrial con res- 
trictasas han puesto de manifesto que 
en las areas en donde coexisten ambas 
formas de Mytilus, estas hibridan en un 
grado variable que depende de las loca- 
lidades (Skibinski, Ahmad y Beard- 
more, 1978; Skibinski y Beardmore, 
1979; Gosling y Wilkins, 1981; Ski¬ 
binski, 1983, 1985; Skibinski, Beard- 
more y Cross, 1983; Edwards y Ski¬ 
binski, 1987; Gosling y McGrath, 1990; 
Gardner y Skibinski, 1991; Coustau et 
al., 1991; McDonald, Seed y Koehn, 
1991; SanjuAn, 1992; ver para una revi¬ 
sion, Gardner, 1992 y Gosling, 1992b). 

Distintos trabajos utilizando genes 
alozfmicos parcialmente diagnostics 
entre Mytilus edulis y M. galloprovmcialis 
han mostrado que en las areas de las Is¬ 
las Britanicas en las que coexisten am- 
bos Mytilus, la composicion genetica de 
estas poblaciones mixtas-hibridas de¬ 


pende de diversos factores ecologies. 
Los alelos caracterfsticos de M. gallopro¬ 
vincialis abundan mas que los de M. edu¬ 
lis en las zonas expuestas al oleaje, 
mientras que, en las zonas resguar- 
dadas, los alelos de M. edulis son los 
predominantes (Gosling y Wilkins, 
1977; 1981; Skibinski et al., 1983). No 
obstante, Skibinski y Roderick (1991) 
mencionan no haber encontrado eviden¬ 
ces claras de este efecto entre zonas ex¬ 
puestas y resguardadas en las localida- 
des que investigan. Varias hipotesis se 
han planteado para tratar de explicar la 
diferencia en frecuencias alelicas entre 
ambientes expuestos y resguardados. 
Skibinski et al. (1983) apuntan hacia 
una mortalidad selectiva de M. edulis 
frente a M. galloprovincialis en el perfodo 
de postfijacion de los mejillones. Con re- 
lacion a esto, Skibinski (1983), Gardner 
y Skibinski (1988) y Skibinski y Rode¬ 
rick (1991) han mostrado que los indivi- 
duos de mayor longitud presentan una 
mayor frecuencia de alelos caracteristi- 
cos de M. galloprovincialis, sugiriendo 
que existe una mortalidad diferencial 
entre ambos Mytilus. No obstante. Gos¬ 
ling y McGrath (1990) no encuentran 
una correlation entre las frecuencias tf- 
picas de M. galloprovincialis para los loci 
Odh y Est-D y la longitud de la concha 
(ver Gosling, 1992b). Por otro lado, 
Gardner y Skibinski (1991) y Willis y 
Skibinski (1992) han puesto de mani- 
fiesto que los individuos de M. gallopro¬ 
vincialis se sujetan con mayor fuerza al 
sustrato que los individuos de M. edulis, 
lo que podrfa explicar la menor mortali- 
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dad de M. galloprovincialis respecto a M. 
edulis en la zona expuesta frente a la 
zona resguardada y, por tanto, una ma¬ 
yor presencia de alelos caracteristicos de 
M. galloprovincialis. Tambien se ha en- 
contrado que los alelos enzimaticos dis- 
tintivos de M. galloprovincialis abundan 
mas en las poblaciones de mejillones fi- 
jados en los niveles altos que en los ba- 
jos en las zonas expuestas de Gran Bre- 
tana (Skibinski, 1983) y quizas tambien 
en las resguardadas (Skibinski y Rode¬ 
rick, 1991). En las costas expuestas de 
Irlanda, Gosling y McGrath (1990) han 
puesto de relieve tambien en el locus 
Odh, pero no en Est-D, una mayor abun- 
dancia de alelos propios de M. gallopro¬ 
vincialis en el nivel alto frente al nivel 
bajo de fijacion. 

El objetivo de este trabajo es caracte- 
rizar geneticamente las poblaciones de 
mejillones procedentes de una localidad 
de la zona hfbrida del sur de la costa atlan- 
tica europea, posiblemente la zona de con- 
tacto entre M. edulis y M. galloprovincialis 
mas antigua segun la hipotesis de Bar- 
sotti y Meluzzi (1968) y a la luz de los 
recientes resultados de SanjuAn et al. 
(1990), Coustau et AL. (1991) y SanjuAn 
(1992). Las muestras proceden de un area 
expuesta al oleaje y otra proxima mas res¬ 
guardada y en cada area, de un nivel u 
horizonte de fijacion alto y de otro bajo, 
con el fin de poner de relieve la inciden- 
cia y el grado de universalidad del efecto 
del microhabitat sobre la composition 
genetica de las poblaciones mixtas e hfbri- 
das de M. edulis y M. galloprovmcialis. 

MATERIAL Y METODOS 

Se estudiaron cinco poblaciones de 
mejillones de la zona hfbrida de Mytilus 
edulis y Mytilus galloprovincialis en la 
costa interna del Golfo de Vizcaya (Cap- 
breton) (Fig. 1). Se recolectaron indivf- 
duos procedentes de un area expuesta 
(Faro de Capbreton) y de otra res¬ 
guardada del oleaje (Canal de Capbre¬ 
ton) separadas unos 200 metros. En 
ambas areas se tomaron muestras al azar 
del nivel u horizonte alto y bajo (Tabla I). 


Con fines comparativos se estudiaron 
tambien dos poblaciones control de M. 
galloprovincialis y una de M. edulis proce¬ 
dentes de la Ria de Vigo (NO de Espana) 
y de Holanda respectivamente. 

Las muestras se recolectaron desde 
abril a julio de 1992, excepto la de Cap¬ 
breton Faro alto (CFA) y la de Holanda 
(ME), que lo fueron en Marzo de 1988 y 
1990 respectivamente. Los ejemplares se 
trasladaron vivos al laboratorio y se con- 
gelaron fntegros a -70 °C hasta su utiliza¬ 
tion. Los mejillones de Rande (GA2) se 
diseccionaron y se les extrajo la glandula 
digestiva que, dividida en dos fragmen- 
tos, sirvio como patron en todo el proceso. 
Se aplico la tecnica de electroforesis hori¬ 
zontal en gel de almidon siguiendo las 
pautas de Harris y Hopkinson (1976) y 
Pasteur, Pasteur, Bonhomme, Catalan 
y Britton-Davidian (1987), y se utilizo la 
glandula digestiva como fuente proteica. 
Se investigaron cuatro loci enzimaticos 
diagnostics, efectivos en la identification 
de las dos formas de mejillon (Ahmad y 
Beardmore, 1976; Skibinski et al., 1978, 
1983; Grant y Cherry, 1985; Varvio, 
Koehn y Vainola, 1988; Beamount, Seed 
y Garcia-Marttnez, 1989; McDonald et 
al., 1991): esterasa-D (Est-D), leucin amino- 
peptidasa-1 (Eap-1), manosa fosfato isomerasa 
(Mpi) y octopi'n deshidrogenasa (Odh) 
siguiendo los protocolos electroforeticos 
y la terminologfa indicados en SanjuAn 
ET AL. (1990). 

A partir de los genotipos electrofore¬ 
ticos observados se calcularon las fre- 
cuencias de los diferentes alelos de cada 
loci. Para analizar la estructura genetica 
de las poblaciones se han calculado los 
estadfsticos F de Wright (1969) siguien¬ 
do a Nei (1977, 1986) y con las correccio- 
nes de Nei y Chesser (1983) que los ha- 
cen independientes del numero de 
muestras o subpoblaciones. La estima 
de F en cada poblacion se realizo utili- 
zando el estadfstico f desarrollado por 
Robertson y Hill (1984). 

Con objeto de sintetizar la information 
obtenida y de permitir un reconocimiento 
mas facil de las relaciones existentes en¬ 
tre las diferentes poblaciones, se han uti- 
lizado las tecnicas estadfsticas multiva- 
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Figura 1. Distribution geografica de Mytilus edulis Linne, 1758 y de M. galloprovincialis 
Lamarck, 1819 en la costa suroccidental europea, indicando la zona de coexistencia e hibrida- 
cion de ambas formas de Mytilus (segun SanjuAn, 1992). Se indican las localidades donde se 
recolectaron las muestras de mejillones. 

Figure L Geographical distribution of Mytilus edulis Linne, 1758 and M. galloprovincialis 
Lamarck, 1819 in the southwestern European coast, showing the coexistance and hybridization 
zone (after SanjuAn, 1992). Localities where mussel sampleswere taken are shown. 


riable (Sneath y Sokal, 1973; Dunn y 
Everitt, 1982). Se ha calculado la distan- 
cia genetica de Nei (1972) entre cada par 
de poblaciones y a partir de la matriz de 
distancias se ha realizado un analisis de 
agrupamientos (UPGMA) y otro de orde- 
nacion (escalado multidimensional). En 
la representation de esta ultima se ha so- 
breimpuesto el arbol expandido mfnimo 
con el fin de visualizar posibles distor- 
siones locales. Para los calculos de los di- 


ferentes analisis se han utilizado los pro- 
gramas informaticos GENET2 (Quesada, 
SanjuAn y SanjuAn, 1992) y NTSYS 
(Rohlf, 1990). 

RESULTADOS 

Las poblaciones analizadas con su 
grado de exposition, nivel de fijacion y 
su situation geografica, se muestran en 
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Tabla I. Poblaciones de mejillones analizadas de Mytilus galloprovincialis (GA), M. edulis 
(ME) y de la zona hfbrida de ambas formas (C), con indicacion del grado de exposicion al ole- 
aje y del nivel del horizonte de recoleccion de las muestras. Cod es un identificador de las 
poblaciones y N es el tamano maximo de la muestra. 

Table I. Studied populations of Mytilus galloprovincialis (GA), M. edulis (ME) and the hybrid 
zone (C), showing the degree of wave action and the collection zones. Cod is a population code 
and N the maximum sample size. 


Poblacion 

Cod 

Localidad 

N 

Exposicion y 
nivel del horizonte 

Holanda 2B 

ME 

Holanda, Cultivo 

34 

Protegida, Bajo 

Capbreton 3A 

C3A 

Capbreton, Canal 

132 

Protegida, Alto 

Capbreton 3BC 

C3BC 

Capbreton, Canal 

133 

Protegida, Bajo 

Capbreton 4A 

C4A 

Capbreton, Faro 

123 

Expuesto, Alto 

Capbreton 4B 

C4B 

Capbreton, Faro 

133 

Expuesto, Bajo 

Capbreton FA 

CFA 

Capbreton, Faro 

133 

Expuesto, Bajo 

Patos B 

GA1 

Patos (Ria de Vigo) 

126 

Protegida, Bajo 

Rande B 

GA2 

Rande (Ria de Vigo) 

133 

Protegida, Bajo 


la Tabla I y en la Figura 1. Las frecuen- 
cias alelicas se presentan en la Tabla II 
agrupando, para cada locus enzimatico, 
los alelos caracterfsticos de Mytilus edulis 
en el alelo compuesto «E», los tfpicos de 
M. galloprovincialis en el «G» y los restan- 
tes en «otros» (vease SanjuAn, 1992). Las 
poblaciones de la zona hfbrida presentan 
frecuencias alelicas de «E» y «G» inter- 
medias entre las poblaciones control de 
M. edulis (ME) y de M. galloprovincialis 
(GA). Las poblaciones de la zona ex- 
puesta (C4A, C4B y CFA) presentan para 
los cuatro loci enzimaticos una mayor 
frecuencia de los alelos G, tfpicos de M. 
galloprovincialis (desde 0,714 a 0,975) que 
las poblaciones de la zona resguardada 
(0,192 a 0,679). Cabe senalar que para to- 
dos los casos, los alelos E, tfpicos de M. 
edulis, se comportan de manera comple- 
mentaria a los de los alelos G. Por otro 
lado, en la zona resguardada, la pobla- 
cion del nivel bajo (C3BC), muestra unas 
frecuencias de los alelos G que oscilan 
desde 0,192 a 0,252, mientras que las del 
nivel alto (C3A) son intermedias entre 
los casos anteriores (0,357 a 0,679). Se de¬ 
tecta asf una significativamente mayor 
frecuencia de los alelos G en el nivel alto 
frente al nivel bajo en la zona res¬ 


guardada (diferencias significativas al 
0,001 en las X 2 de homogeneidad para los 
cuatro loci). En las poblaciones de la zona 
expuesta muestreada en la misma epoca 
(C4), no se detectan diferencias significa¬ 
tivas de las frecuencias alelicas mediante 
la X 2 de homogeneidad entre la muestra 
del nivel alto (C4A) y la del nivel bajo 
(C4B), aun cuando si se puede observar 
una ligera pero consistente mayor pro- 
porcion de los alelos G para los cuatro 
loci en el nivel alto (C4A) que en el bajo 
(C4B) (el test no parametrico de Wilco- 
xon da diferencias significativas al 5%). 

En la Tabla III figuran las estimas de 
las desviaciones de las proporciones del 
equilibrio Hardy-Weinberg (HWE) por 
locus, calculadas mediante el estadfstico 
f de Robertson y Hill (1984) con su 
error estandar, asf como el test X 2 de 
bondad de ajuste al HWE. El examen de 
estas estimas pone de relieve que la 
mayorfa de las poblaciones de la zona 
hfbrida (C) presentan un deficit de hete- 
rocigotos (valores positivos de f), indica- 
tivos de una posible mezcla de poblacio¬ 
nes con frecuencias alelicas distintas 
(efecto Wahlund). La Tabla IV muestra 
las estimas de los indices de fijacion de 
Wright (1969) segun Nei (1977, 1986) y 
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Tabla II. Frecuencias de los alelos tfpicos de Mytilus edulis (E) y M. galloprovincialis (G) de 
los loci diagnostico entre ambas formas de Mytilus en las poblaciones analizadas. N es el tama- 
no de la muestra. 

Table II. Frequencies of the typical alleles o/Mytilus edulis (E) and M. galloprovincialis (G) at 
the diagnostic loci between both forms of Mytilus in the studied populations. N is the sample 
size. 


Locus 

ME 

C3A 

C3BC 

Poblacion 

C4A C4B 

CFA 

GA1 

GA2 

EstD 

G. 

0,074 

0,581 

0,192 

0,781 

0,729 

0,803 

0,894 

0,902 

E 

0,897 

0,378 

0,771 

0,174 

0,244 

0,159 

0,043 

0,034 

otros 

0,030 

0,040 

0,038 

0,045 

0,028 

0,039 

0,063 

0,064 

N 

34 

123 

133 

121 

129 

132 

126 

133 

Lap 7 

G 

0,030 

0,608 

0,252 

0,820 

0,759 

0,763 

0,891 

0,866 

E 

0,956 

0,378 

0,748 

0,149 

0,220 

0,236 

0,086 

0,101 

otros 

0,015 

0,013 

0,000 

0,030 

0,020 

0,000 

0,024 

0,032 

N 

34 

115 

119 

117 

127 

93 

64 

123 

Mpi 

G 

0,044 

0,679 

0,252 

0,858 

0,792 

0,975 

0,938 

0,970 

E 

0,941 

0,321 

0,748 

0,142 

0,208 

0,025 

0,049 

0,017 

otros 

0,015 

0,000 

0,000 

0,000 

0,000 

0,000 

0,012 

0,013 

N 

34 

120 

129 

123 

120 

20 

81 

117 

Odh 

G 

0,031 

0,537 

0,245 

0,747 

0,714 

0,737 

0,851* 

0,818 

E 

0,938 

0,439 

0,748 

0,222 

0,252 

0,256 

0,117 

0,174 

otros 

0,031 

0,024 

0,008 

0,030 

0,035 

0,006 

0,032 

0,008 

N 

32 

123 

131 

117 

133 

78 

94 

132 


Nei y Chesser (1983) (F'is, F'it y F'st) 
para las poblaciones de la zona hibrida 
(C). Los valores relativamente altos de 
F'st (desde 0,0975 a 0,3483) y significati- 
vos al 0,1% (X 2 de homogeneidad) para 
los loci diagnostico, son indicativos de 
un grado de diferenciacion importante 
en una distancia de menos de 200 m. 
Estos valores de F'st elevados y las fre¬ 
cuencias alelicas intermedias en dife- 
rente grado entre las que presentan las 
poblaciones puras de M. edulis y M. 
galloprovincialis ponen de manifiesto la 
mezcla de los dos acervos geneticos en 
una proportion que varia con el micro¬ 
habitat. 

A partir de la matriz de distancias 
geneticas (Nei, 1972) se ha construido un 
dendrograma UPGMA (Fig. 2) donde se 


observan dos grupos bien diferenciados: 
uno formado por la poblacion control de 
M. edulis (ME) y la poblacion de la zona 
resguardada baja (C3BC) y otro con las 
dos poblaciones control de M. gallopro¬ 
vincialis (GA_) y el resto de poblaciones 
de la zona hibrida. En este segundo 
grupo se separan las poblaciones control 
de M. galloprovincialis (GA) de las de la 
zona hibrida del area expuesta (C), 
situandose en posicion aparte la pobla¬ 
cion de la zona resguardada alta (C3A). 
Esta relation se clarifica en la disposition 
espacial que se puede observar en el 
escalado multidimensional con sobreim- 
presion del arbol expandido rrunimo que 
se muestra en la Figura 3, en la que se 
observa que las poblaciones de la zona 
hibrida presentan una posicion interme- 
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Tabla III. Estimas de las desviaciones de las proporciones de equilibrio Hardy-Weinberg (HWE) 
para los loci diagnostico Est-D , Lap-1, Mpi y Odh en las poblaciones de Mytilus edulis (ME), 
M. galloprovincialis (GA) y de la zona hibrida de ambas formas (C). f es la estima de 
Robertson y Hill (1984) con su error estandar (f±SE). % 2 es la prueba de ji-cuadrado de ajus- 
te a las proporciones esperadas en HWE. 

Table III. Estimation of the deviations from the equilibrium Hardy-Weinberg (HWE) proportions 
at the diagnostic loci Est-D, Lap-1, Mpi and Odh in the populations o/Mytilus edulis (ME), M. 
galloprovincialis (GA) and the hybrid zone (C).fis the Robertson and Hill's (1984) estima¬ 
te, together with its standard error (f±SE). x 2 is the chi-square test of fit to the HWE expected 
proportions. 


Locus 

ME 

C3A 

C3BC 

Poblaciones 

C4A C4B 

CFA 

GA1 

GA2 

Est-D 

f 

-0,071 

0,277** 0,234** 

0,140 

0,423*** 

-0,004 

-0,048 

-0,047 

SE 

±0,172 

±0,090 

±0,087 

±0,091 

±0,088 

±0,087 

±0,063 

±0,061 

X 2 

0,45 

12,84* 

14,30 ** 

4,03 

25,17*** 

1,28 

1,78 

1,56 

Lap-1 

f 

0,242 

0,123* 

0,243*** 

0,051 

0,019 

0,159* 

0,512* 

** 0,194** 

SE 

±0,172 

±0,054 

±0,053 

±0,065 

±0,051 

±0,073 

±0,088 

±0,064 

X 2 

2,46 

20,08 ** 

30,30 *** 

8,96 

10,08 

8,51 

38,49*** 15,98 * 

Mpi 

f 

-0,001 

0,144 

0,200* 

0,185* 

-0,006 

0,000 

-0,048 

0,328*** 

SE 

±0,172 

±0,091 

±0,088 

±0,090 

±0,091 

±0,224 

±0,111 

±0,065 

X 2 

- 

2,34 

4,94* 

3,97* 

0,01 

0,01 

0,35 

8,92 ** 

Odh 

f 

0,507*** 

0,142* 

0,245*** 

0,178** 

0,213*** 

0,153 

0,091 

0,107 

SE 

±0,125 

±0,064 

±0,062 

±0,065 

±0,061 

±0,080 

±0,073 

±0,062 

X 2 

co 

o 

CO 

12,25 

16,88*** 

25,13** 

16,78* 

5,38 

4,58 

4,89 

* p < 0,05 

** p < 0,01 

— 

p < 0,001 (Prueba de significacion para f=0, Robertson y 

Hill, 1984). 


Tabla IV. Estimas de los indices de fijacion de Wright (F’is F’it y F’st) (segun Nei y Chesser, 
1983) para las poblaciones de la zona hfbrida de Mytilus galloprovincialis y M. edulis. yf es el 
valor de la prueba ji-cuadrado de homogeneidad. 

Table IV. Estimates of the Wright's fixation indices (F'js F'jt y F'st) (after Nei and Chesser, 
1983) in the populations of the hybrid zone o/Mytilus galloprovincialis and M. edulis. x 2 is the 
homogeneity chi-square value. 

locus 

Fis 

F’it 

F’st 

X 2 

Est-D 

0,1959 

0,3923 

0,2443 

326,9*** 

Lap-1 

0,1194 

0,2053 

0,0975 

252,7*** 

Mpi 

0,1183 

0,4254 

0,3483 

267,8*** 

Odh 

0,1670 

0,2587 

0,1101 

210,0*** 

* p < 0,05 

** p < 0,01 *** p < 0,001 (Nivel de significacion para la prueba de % 2 

de homogeneidad). 
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Distancia genetica 



GA1 

GA2 

CFA 

C4A 

C4B 

C3A 

C3BC 

ME 


Figura 2. Dendrograma UPGMA de las poblaciones de mejillon de la zona hibrida (C) y de las 
poblaciones control de Mytilus galloprovincialis (GA) y de M. edulis (ME) utilizando la matriz 
de distancias geneticas (ver Nei, 1972) obtenida de las frecuencias alelicas de los loci Est-D , 
Lap-1 , Mpi y Odh. El coeficiente cofenetico es de r= 0,829. 

Figure 2. UPGMA dendrogram of mussel populations from the hybrid zone (C) and Mytilus 
galloprovincialis (GA) and M. edulis (ME) control populations, using the genetic distance 
matrix (after Nei, 1972) obtained from the allele frequencies of Est-D, Lap-1, Mpi and Odh loci. 
Cophenetic coefficient is r= 0.829. 



Figura 3. Escalado multidimensional no-metrico, con sobreimpresion del arbol expandido mmi- 
mo, de poblaciones de la zona hibrida (C) y de poblaciones control de M. galloprovincialis 
(GA) y de M. edulis (ME) utilizando la matriz de distancias geneticas (ver Nei, 1972) obtenida 
a partir de las frecuencias alelicas de los loci Est-D , Lap-1, Mpi y Odh. El “stress” final es S= 
0,0072, lo que supone un ajuste excelente o perfecto. 

Figure 3. Non-metric multidimensional scaling, with the minimum spanning tree superimposed , 
of the populations from the hybrid zone (C) and Mytilus galloprovincialis (GA) and M. edulis 
(ME) control populations using the genetic distance matrix (after Nei, 1972) obtained from the 
allele frequencies o/Est-D, Lap-1, Mpi and Odh loci. Final “stress” is S= 0.0072, which means 
and excellent or perfect fit. 
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dia entre las extremas tipicas M. edulis 
(ME) y M. gallopro-vincialis (GA), ubican- 
dose la poblacion de la zona res- 
guardada baja (C3BC) proxima a la de 
M. edulis (ME), las de la zona expuesta se 
agrupan hacia las de control de M. gallo¬ 
provincialis (GA) y en posicion interme¬ 
dia se situa la poblacion de la zona res- 
guardada alta (C3A). 

DISCUSION 

El analisis de los cuatro loci enzimati- 
cos diagnostico ( Est-D, Lap-1, Mpi y Odh) 
en las cinco poblaciones de la zona 
hfbrida muestra importantes diferencias 
en las frecuencias alelicas de cada pobla¬ 
cion estudiada, segun el habitat de 
donde proceden. Las poblaciones del 
area expuesta (C4A, C4B y CFA) presen- 
tan las frecuencias alelicas para los 
cuatro loci semejantes a las tipicas de M. 
galloprovincialis (Tabla II), agrupandose 
asi con las dos poblaciones control de M. 
galloprovincialis (GA) tanto en el dendro- 
grama UPGMA (Fig. 2) como en el esca- 
lado multidimensional (Fig. 3). En las 
poblaciones de la zona resguardada, la 
del nivel bajo (C3BC) presenta una 
elevada y significativa frecuencia de los 
alelos tipicos de M. edtdis para los cuatro 
loci estudiados y la poblacion del nivel 
alto (C3A) mantiene unas frecuencias 
alelicas claramente intermedias entre las 
de las poblaciones control de M. gallopro¬ 
vincialis (GA) y las de M. edulis (ME) 
(Tabla II). El dendrograma UPGMA y el 
escalado multidimensional agrupan a la 
poblacion resguardada del nivel bajo 
(C3BC) con la control de M. edulis, ubi- 
candose en posicion intermedia la del 
nivel alto (C3A) (Figs. 3 y 4). Estos datos 
ponen de manifiesto la mayor propor- 
cion de los alelos tipicos de M. gallopro¬ 
vincialis (G) que los de M. edulis (E) para 
los loci Est-D, Lap-1, Mpi y Odh en la zona 
expuesta frente a la zona resguardada. 
Se corrobora asi este efecto encontrado 
por otros autores en las Islas Britanicas 
para Gpi, Lap-2 y Pgm (Gosling y 
Wilkins, 1977 y 1981) y para Est-D (Ski- 
binski ET al., 1983), aun cuando Ski- 
binski y Roderick (1991) manifiesten no 


haber encontrado pruebas claras de 
dicho efecto. Para tratar de explicar este 
comportamiento diferencial entre pobla¬ 
ciones de zonas expuestas y resguar- 
dadas, Skibinski et al. (1983) han suge- 
rido una mortalidad selectiva menor de 
los juveniles de M. galloprovincialis frente 
a los de M. edulis . Con relacion a esto, 
Skibinski (1983), Gardner y Skibinski 
(1988) y Skibinski y Roderick (1991) han 
encontrado que la frecuencia de los 
alelos «G» se incrementa con el tamano 
de las conchas, indicando que esta corre- 
lacion bien pudiera deberse a una morta¬ 
lidad diferencial mayor de los indivi- 
duos juveniles de M. edulis, debido posi- 
blemente a la accion de las olas o al 
efecto erosivo de la arena suspendida en 
el agua. En lfnea con esta argumenta- 
cion, varios estudios han mostrado que, 
en el mismo habitat, los individuos M. 
edulis se separan con mayor facilidad del 
sustrato que los M. galloprovincialis 
(Gardner y Skibinski, 1991; Willis y 
Skibinski, 1992). Esta distinta capacidad 
de sujecion a la roca podria actuar como 
una diferencia adaptativa causante de 
una mayor mortalidad de los M. edulis 
juveniles frente a los de M. galloprovincia¬ 
lis en las zonas expuestas. Asi se podria 
explicar la mayor frecuencia de alelos 
tipicos de M. galloprovincialis que de M. 
edulis en estas zonas. No obstante, 
alguna ventaja adaptativa debe poseer 
M. edulis frente a M. galloprovincialis para 
contrarrestar los efectos de su mayor 
mortalidad y su consistentemente menor 
fertilidad (Gardner y Skibinski, 1990), y 
poder mantener su identidad genetica de 
manera aparentemente estable en el dis- 
currir de los anos en zonas hfbridas (en 
Rock durante 30 anos, segun Beaumont, 
Seed y Garcia-MartInez, 1989). 

Por otro lado, los resultados de este 
trabajo muestran diferencias en las fre¬ 
cuencias alelicas entre muestras proce- 
dentes de distintos niveles de fijacion. 
En la zona resguardada, los alelos de los 
cuatro loci Est-D G , Lap-1 G , Mpf y Odh G , 
abundan mas en el nivel alto (entre 50 y 
70%) que en el nivel bajo, en donde pre- 
dominan significativamente los alelos E, 
(alrededor de un 75%). Sin embargo, en 
la zona expuesta, no se aprecian diferen- 


109 





IBERUS , 14 (2), 1996 


cias significativas en las frecuencias ale- 
licas entre ambos niveles de fijacion, aun 
cuando se observa la tendencia de una 
mayor frecuencia de los alelos G que los 
E (sobre un 5% mas en cada locus) en el 
nivel superior frente al inferior. Estos re- 
sultados, si bien confirman la idea gene¬ 
ral de una mayor frecuencia de alelos G 
en los niveles superiores frente a los in¬ 
feriors, por otro lado contrastan relati- 
vamente con lo encontrado en las zonas 
hibridas de las Islas Britanicas. En zonas 
expuestas, Skibinski (1983) ha encon¬ 
trado en Gran Bretana una mayor fre- 
cuecia de los alelos Est-D G f Lap-1 G y 
Odh G en los niveles altos frente a los ba- 
jos, mientras que Gosling y McGrath 
(1990), en la costa Oeste de Irlanda, tan 
solo han detectado diferencias significa¬ 
tivas en el locus Odh pero no en Est-D. 
No obstante, el trabajo de Skibinski y 
Roderick (1991) parece sugerir que este 
efecto en Est-D se ha detectado tanto en 
zonas resguardadas como en expuestas 
de Gran Bretana. Por otro lado, la ante- 
riormente senalada correlacion entre la 
frecuencia de los alelos caracteristicos 
de M. galloprovincialis y la edad de los 
individuos (Skibinski, 1983; Skibinski y 
Roderick, 1991), junto con la interac- 
cion de los distintos factores ambienta- 
les condiciona la comprension de la 
composicion genetica y la dinamica de 
la poblaciones mixtas de Mytilus. Asi, 
Gosling y McGrath (1990) han encon¬ 
trado, en dos localidades expuestas, di¬ 
ferencias en las frecuencias alelicas de 
mejillones del nivel medio y del bajo 
mostrando respectivamente una fre¬ 
cuencia del alelo Odh G de 0,612 y 0,544 
en Ballynahown y de 0,549 y 0,414 en 
Carraroe. En el presente trabajo, en la 
zona expuesta de Capbreton, la frecuen¬ 
cia de Odh G es de 0,747 en el nivel alto y 
de 0,714 en el bajo. La comparacion de 
los datos indica que en Capbreton hay 
un porcentaje considerablemente mayor 
de los alelos Odh G que en Ballynahown 
o Carraroe; esta diferencia podria ser in- 
dicativa de que el agente selectivo (por 
ejemplo, la exposicion al oleaje actuando 
sobre la diferente fuerza de fijacion al 
sustrato de M. edulis y M. galloprovincia¬ 


lis) esta incidiendo mas intensamente en 
Capbreton que en las dos poblaciones ir- 
landesas, si se considera la correlacion 
encontrada entre la intensidad de expo¬ 
sicion al oleaje y la frecuencia del alelo 
Est-D G ( =Est-D 8 ; Skibinski et al ., 1983). 
Ademas, bien podria pensarse que la in¬ 
tensidad del agente selectivo puede lie- 
gar a ser tan elevada que otros factores 
(como los que contribuyen a establecer 
las diferencias de habitat alto /bajo) po- 
drfan disminuir relativamente en impor- 
tancia, de modo que el primer factor es- 
tarfa homogeneizando decisivamente el 
habitat de la zona expuesta anulando las 
diferencias geneticas correlativas entre 
las poblaciones mixtas de Mytilus de ni¬ 
veles altos y bajos. La utilizacion de una 
escala de exposicion al oleaje podria 
ayudar a comprender en un futuro la in- 
teraccion de los distintos factores selecti- 
vos. Por tanto, hasta el momento, no pa¬ 
rece existir una clara concordancia en la 
constatacion y delimitacion de dicho fe- 
nomeno, precisando seguramente de 
mas analisis que tengan en cuenta la 
edad de los individuos y/o larvas y la 
amplitud e intensidad de los factores de 
los distintos ambientes objeto de estu- 
dio. 

No obstante, varias hipotesis se es- 
tan barajando para dar cuenta de esta 
microdiferenciacion. Como explicaron 
Gosling y McGrath (1990), las diferen¬ 
cias geneticas entre los mejillones de 
distintos niveles de fijacion pueden de- 
berse a cohortes geneticamente diferen- 
tes de las larvas que se fijan preferente- 
mente en distintos niveles de la costa o 
bien a que las larvas que Began a la 
costa sean geneticamente homogeneas y 
una vez fijadas diverjan geneticamente 
por la presion selectiva de distintos 
agentes como por ejemplo, por el dife¬ 
rente grado de exposicion al aire segun 
el nivel de fijacion. Skibinski, (1983) y 
Skibinski y Roderick (1991) se inclinan 
por esta ultima hipotesis sosteniendo 
que M. galloprovincialis esta mejor adap- 
tado que M. edulis a sujetarse al sustrato 
y a explotar mas favorablemente el es- 
pacio disponible por razon de su pecu¬ 
liar morfologfa. Sm embargo, las dife- 
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rencias genicas en Est-D y Lap-1 entre 
cohortes de larvas que se fijaron en 
cuerdas en momentos distintos en el pe- 
riodo de 2 semanas (Gosling y Wilkins, 
1985), y las diferencias genicas en Est-D 
y Odh encontradas entre mejillones fija- 
dos en distintos niveles (Gosling y Mc¬ 
Grath, 1990) apuntan hacia diferencias 
geneticas originadas ya en la fase larva- 
ria o bien en los primeros estadios des¬ 
pues de la fijacion (hasta un mes). Por 
otro lado, en el caso de diferencias en es- 
tadfos tan juveniles, no parece aceptable 
pensar, por razon del pequeno tamano 
de ambos Mytilus, en una ventaja adap- 
tativa morfologica de M. galloprovincialis 
frente a M. edulis con respecto al espa- 
cio; ademas, en los habitats en los que 
abundan los alelos G es donde tiene lu- 
gar una mayor convergencia morfolo¬ 
gica de la forma fenotipicamente mas 
plastica M. galloprovincialis hacia la mas 
homogenea M. edulis (Seed, 1978; Gos¬ 
ling, 1984). 

Sean cuales sean los factores ambien- 
tales que conforman las diferencias gene¬ 
ticas de los ambientes expuesto / resguar- 
dado y nivel alto/bajo, estos parecen in- 
teraccionar entre sf de manera que la 
composition genetica de las poblaciones 
mixtas de Mytilus no resulta facilmente 
previsible. Sin embargo, parece claro que 
estos diferentes habitats deberan tenerse 
en cuenta a la hora de tratar de entender 
los distintos aspectos de la dinamica de 
las zonas hibridas de ambos Mytilus, 
como por ejemplo la ausencia de gra- 
dientes genicos claros en la zona de con- 
tacto en la costa atlantica continental 
(COUSTAU ET AL., 1991). Valga la conside¬ 
ration de que precisamente el area inves- 
tigada en el presente trabajo bien pu- 
diera ser la zona de contacto de M. edulis 
y M. galloprovincialis mas antigua segun 
la hipotesis de Barsotti y Meluzzi 
(1968) en la que se supone que la presen- 
cia de M. galloprovincialis en las costas 
noroccidentales de Europa se atribuye a 
una migration desde el Mar Mediterra- 
neo. Recientemente SanjuAn et al., 
(1986, 1990), Coustau et al., (1991) y 
SanjuAn (1992) han establecido inequi- 
vocamente una distribution continua de 


M. galloprovincialis desde el Mar Medite- 
rraneo hasta el Canal de la Mancha y las 
Islas Britanicas, haciendo mas plausible 
la hipotesis anterior (ver Gosling, 
1992a). Por tanto, cabria esperar una am- 
plia homogeneizacion genetica en ese 
area mas antigua, bien por una impor- 
tante introgresion de ambos Mytilus, o 
bien por eliminacion de M. edulis por se¬ 
lection direccional. Sin embargo, en Cap- 
breton se encuentra un porcentaje im- 
portante de alelos G en la zona expuesta 
(70-85%), mientras que en la zona res- 
guardada del nivel bajo tan solo alcanza 
el 25% y los alelos E oscilan sobre el 75%. 
Si se admite en esta area una mayor mor- 
talidad y una menor fertilidad de M. edu¬ 
lis respecto de M. galloprovincialis (como 
lo descrito en otras areas por Gardner y 
Skibinski, 1990) no resultan faciles de ex- 
plicar las diferencias genicas entre zonas 
tan proximas, maximo si se considera 
que el sentido hacia el norte de las co- 
rrientes marinas superficiales predomi- 
nantes redundaria en un aporte masivo 
de larvas de M. galloprovincialis. De 
nuevo se suscita el problema de que al- 
guna ventaja adaptativa deben de pre- 
sentar los individuos de M. edulis para 
permanecer en una proporcion no 
desdenable en un area que, suponiendo 
actuaran los factores hasta ahora conoci- 
dos, resultana en un proceso selectivo 
unidireccional con eliminacion de M. 
edulis y de sus alelos caracteristicos. jY 
esto no parece ser lo que sucede! 
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